System size dependence of cluster properties from two-particle angular correlations in 
Cu+Cu and Au+Au collisions at ^/s^ = 200 GeV 
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We present results on two-particle angular correlations in Cu-|-Cu and Au-f Au collisions at a center 
of mass energy per nucleon pair of 200 GeV over a broad range of pseudorapidity (rj) and azimuthal 
angle {(j)) as a function of collision centrality. The PHOBOS detector at RHIC has a uniquely- 
large angular coverage for inclusive charged particles, which allows for the study of correlations 
on both long- and short-range scales. A complex two-dimensional correlation structure in Arj and 
Atf) emerges, which is interpreted in the context of a cluster model. The effective cluster size and 
decay width are extracted from the two-particle pseudorapidity correlation functions. The effective 
cluster size found in semi-central Cu-|-Cu and Au-|-Au collisions is comparable to that found in 
proton-proton collisions but a non-trivial decrease of the size with increasing centrality is observed. 
Moreover, a comparison between results from Cu-|-Cu and Au-|-Au collisions shows an interesting 
scaling of the effective cluster size with the measured fraction of total cross section (which is related 
to the ratio of the impact parameter to the nuclear radius, b/2R), suggesting a geometric origin. 
Further analysis for pairs from restricted azimuthal regions shows that the effective cluster size at 
Aff) ~ 180° drops more rapidly toward central collisions than the size at A(j!) ~ 0°. The effect of 
limited rj acceptance on the cluster parameters is also addressed, and a correction is applied to 
present cluster parameters for full rj coverage, leading to much larger effective cluster sizes and 
widths than previously noted in the literature. These results should provide insight into the hot 
and dense medium created in heavy ion collisions. 

PACS numbers: 25.75.-q,25.75.Dw,25.75.Gz 



I. INTRODUCTION 

Multiparticle correlation studies have played an im- 
portant role in exploring the underlying mechanism of 
particle production in high energy hadronic collisions. In 
p-|-p collisions, inclusive two-particle correlations have 
been found to have two components: "intrinsic" two- 
particle correlations as well as an effective "long-range" 
correlation due to event- by- event fluctuations of the 
overall particle multiplicity (see, for example, Ref. 
By considering the two-particle rapidity density at fixed 
multiplicity, the intrinsic correlations between particles 
were isolated and found to be approximately Gaussian, 
with a range of ^ 1 unit in pseudorapidity (ry = 



— ln(tan(6'/2))). These correlations have been conven- 
tionally called "short-range" . Their properties have been 
characterized by the concept of "cluster emission" [2] . 

The simple idea that hadrons are produced in clus- 
ters, rather than individually, has had great success in 
describing many features of particle production . In 
this scenario, hadronization proceeds via "clusters" , high 
mass states (e.g. resonances but not necessarily with all 
well-defined quantum numbers) which decay isotropically 
in their rest frame into the final-state hadrons. An 
independent cluster emission model (ICM) has been 
widely applied to the study of two-particle correlations 
[ll-Q, where clusters are formed before the final-state 
hadrons and are independently emitted according to a 
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dynamically generated distribution in i] and (j). The 
clusters subsequently decay isotropically in their own rest 
frame into the observed final-state hadrons. The ob- 
served correlation strength and extent in phase space can 
be parameterized in terms of the cluster multiplicity, or 
"size" (the average number of particles in a cluster) and 
the decay "width" (which characterizes the separation 
of the particles in pseudorapidity) . However, it should 
be noted that independent cluster emission is only a 
phenomenological approach which provides no insight as 
to the mechanisms by which clusters are formed. Further 
modeling is required to connect these studies to the 
underlying QCD dynamics. 

A measurement of cluster properties from two-particle 
correlations in p+p collisions for particles emitted into 
\r]\ < 3 was performed previously at center of mass 
energies {^/s) of 200 GeV and 410 GeV using the 
PHOBOS detector at the Relativistic Heavy Ion Collider 
(RHIC) The data suggested an effective cluster size 
{K^ff, defined in Sect.|IV| of2.44±0.08 at Vs = 200 GeV, 
increasing with collision energy and event multiplicity. 
The results are consistent with previous measurements 
from ISR and SppS at various energies [H, H, |1] ■ 

In heavy ion collisions at RHIC, it has been predicted 
that the formation of a strongly interacting quark gluon 
plasma could modify cluster properties relative to p-l-p 
collisions [7,]. In order to systematically explore the 
properties of the clusters from p-t-p to A-l-A collisions, 
this paper presents the results on two-particle angular 
correlations in Cu-|-Cu and Au-f Au collisions at center 
of mass energy per nucleon pair (y^s^) of 200 GeV, 
over a very broad acceptance in A77 (= ?7i — '72) and 
(= 01-02). The PHOBOS Octagon detector, 
covering pseudorapidity -3 < 77 < 3 over almost the 
full azimuth, is well suited to measure the correlations 
between particles emitted from clusters. The two- 
dimensional (2-D) correlation functions in Cu+Cu and 
Au-I-Au collisions, as well as the extracted effective 
cluster size and width from one-dimensional (1-D) Arj 
correlation functions, are presented as a function of 
system size. By separating particle pairs into "near- 
side" (0° < A0 < 90°) and "away-side" (90° < A0 < 
180°), more detailed information on the properties of the 
clusters is obtained. Furthermore, extrapolating limited 
acceptance (It^I < 3) to full phase space, cluster properties 
unbiased by detector acceptance have been estimated. 
This comprehensive analysis of cluster properties in p-l-p 
and A-l-A collisions should provide useful information for 
understanding the hadronization stage, but may also give 
insight into physics relevant at much earlier times. 



II. DATA SETS 

The data presented here for Cu-l-Cu and Au-f Au 
collisions at ^s„„ = 200 GeV were collected during 
RHIC Run 4 (2004) and Run 5 (2005) using the large- 
acceptance PHOBOS Octagon silicon array covering 
pseudorapidity —3 < < 3 over almost the full azimuth. 
A full description of the PHOBOS detector can be found 
in Ref. [8]. The primary event trigger used the time 
difference between signals in two sets of 10 Cerenkov 
counters located at 4.4 < \ri\ < 4.9, to select collisions 
that were close to the nominal vertex position Zytx — 
along the beam-axis. About 4 million events each of 
Cu-I-Cu and Au-|-Au colhsions at y/J^ = 200 GeV were 
selected for further analysis by requiring that the primary 
collision vertex falls within | Zytx \ < 6 cm. 

The angular coordinates (rj, (f>) of charged particles are 
measured using the location of the energy deposited in 
the silicon pads of the Octagon. The granularity of the 
Octagon is determined by the sizes of the readout pads 
which are about 11.25° (~ 0.2 radians) in and range 
from 0.006 to 0.05 in 77. Noise and background hits are 
rejected by placing a lower threshold on the deposited 
energy corrected for the path length through the silicon 
after hit merging, assuming that the charged particle 
originated from the primary vertex. Depending on 77, 
merged hits with less than 50-60% of the energy loss 
expected for a minimum ionizing particle are rejected. 
More details of the hit reconstruction procedure can be 
found in Ref. [9j. 

III. ANALYSIS PROCEDURE 

The detailed analysis procedure is described in 
Ref. [6|]. The inclusive two-particle correlation function 
in (A77, A(j)) space is defined as follows: 

^ ((n-l)p"(Ar;,A0)) ^ _ 

The first line of Eq. [T] shows the behavior of the correla- 
tion function, while its mathematically equivalent second 
line corresponds to the analysis procedure described 
below. For a given centrality bin, p^^{Ar], Acf)) (normal- 
ized to unit integral) is the foreground pair distribution 
evaluated event- by- event by taking pairs within the same 
event. It is then weighted by a factor, (n — 1) (where n 
is the total number of hits in each event, while generator 
level tracks are used in MC), and averaged over all the 
events. The event multiplicity normalization factor, (n — 
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FIG. 1: (a) dE/dx distribution for hits registered in the Octagon in the < r) < 0.3 range (open circles) fitted by a sum (solid 
line) of predicted dE/dx distributions for various number of particles hitting a single pad (dashed lines) for the most central 
3% of 200 GeV Au+Au collisions. A MIP is defined as the energy deposited by a single minimum-ionizing particle at normal 
incidence, (b) Estimated average number of particles per silicon pad as a function of rj and dE/dx for the most central 3% of 
200 GeV Au+Au collisions. 



1), is introduced to compensate for the trivial dilution 
effects from uncorrelated particles since the number of 
uncorrelated pairs grows quadratically with n, while the 
number of correlated pairs grows only linearly. Therefore, 
R{Ar], A(j)) is defined in such a way that if a heavy 
ion collision is simply a superposition of individual p+p 
collisions, the same correlation function will be observed 
in both A+A and p+p collisions. 

The mixed-event background distribution 
/9™'™'^(A77, A(j)) (also normalized to unit integral) is 
constructed by randomly selecting two particles from 
two different events with similar vertex (with a bin 
width of 0.2 cm) and centrality (in bins of 5% of the 
total inelastic cross section), representing a product of 
two single particle distributions. Since the background 
is found to be roughly multiplicity independent within 
a centrality bin, the inclusive p'^^^^'^ {Arj , Acj)) is used 
in our calculations of Eq. [T] for each centrality class. 
Thus it can be either inside or outside the event average 
as shown in the two lines of Eq. [TJ R{Ari, Acj)) is 
symmetrized at A?7=0 and A(l)—0 following the same 
convention as in the previous publication j6|. As can 
be inferred from the first line of Eq. [U since both 
p^^{Ar], Acj)) and yO™"'°'^(A?7, A0) are normalized to unit 
integral, R{Ar], Ac/)) will range from positive to negative. 
Therefore, there is no particular significance when 
R{Ar], A(/))=0 in a small region of Ary and A(j), although 
R{Ar], A(j))=0 everywhere would indicate the absence of 
correlations. 

Corrections for secondary effects and incomplete ac- 
ceptance have been applied in a similar way to those in 



p+p collisions resulting in corrections with comparable 
magnitudes 6]. Since the PHOBOS Octagon is a single- 
layer silicon detector, there is no pr, charge or mass infor- 
mation available for the particles. All charged particles 
above a \ow-pT cutoff of about 7 MeV/c at 77=8, and 
35 MeV/c at 77=0 (which is the threshold below which 
a charged particle is stopped by the beryllium beam 
pipe) are included on equal footing. Thus secondary 
effects, such as ^-electrons, 7 conversions and weak 
decays, cannot be all rejected directly. The incomplete 
azimuthal acceptance in some pseudorapidity regions 
naturally suppresses the overall correlation strength, but 
MC simulations show that it does not change the shape 
of the correlation function R{Ari, Acj}). To correct for 
these detector effects in the data, correlation functions 
are calculated for MC events (PYTHIA for p+p collisions 
and HIJING for A+A coUisions) at = 200 GeV both 
at the generator level for true primary charged hadrons, 
i?^.p(A7?, A(f>), and with the full GEANT detector simu- 
lation and reconstruction procedure, R^^'^lArj, A<j)). The 
whole correction procedure can be summarized by the 
following equation; 

Rt'^{Ar,,A^) = Ax [i?da;a(A^, A0)-5(Ar?, A</.)]. (2) 

where the function, 

SiAr,, Acj}) = R^^{A^, Acj>) ~ <p,ec(Ar?, A^) (3) 

represents the correction for secondary effects which 
mainly contribute to a large, narrow spike in R{Ari, Acj>) 
around A7/=0, Acj>=0 (for more details, see Ref. [6]). 
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The acceptance correction factor, A, is estimated to 
be around 1.3 with little dependence on vertex and 
centrality. The total systematic uncertainty for these 
corrections is typically about 5% of the final value of 
i?(A77,A0). 

In addition to these corrections, the high occupancies 
measured in A+A collisions (which are in the range 
from 50-60% at midrapidity for the most central 3% 
of 200 GeV Au+Au collisions) require us to account 
for the high probability of multiple particles hitting a 
single pad. The dE/dx distribution of hits in a very low 
multiplicity environment (e.g. 55-60% peripheral Cu+Cu 
which has an occupancy of about 4% at midrapidity) 
has been measured first in a narrow bin of 0.3 unit 
(same bin size as for i?(A77, A0)), to approximate the 
dE/dx distribution for a single particle hitting a single 
silicon pad, Piijgr)- By convoluting Pi(^) i times, 
the dE/dx distribution of i particles hitting a single 
pad can be predicted, Pi{^) (i=l, 2, 3....). Then, the 
dE/dx distribution in more central data events, p{^), 
is fitted by a sum of Pi{^) with weighting factor wi 
such that the relative contribution of different number 
of particles hitting on a single pad can be estimated, as 
illustrated in Fig. [1^ for the most central 3% of 200 GeV 
Au-|-Au collisions as an example. While the Wi are free 
parameters in these fits, they have been compared to 
a Poisson distribution and found to be in very good 
agreement. The average number of particles per pad can 
thus be calculated as a function of dE/dx, 



dE 



(4) 



Performing the procedure above in each rj and centrality 
bin, the average number of particles per pad as a function 
of dE/dx and rj, N{ri, ^), is derived and shown in 
Fig. [T}d for the most central 3% of 200 GeV Au+Au 
collisions. In calculating the correlation function, each 
hit is assigned a weight N{?], ^) based on its ?/ and 
dE/dx. Using this procedure, the effects of high occu- 
pancy are compensated at the hit level. The systematic 
uncertainty of the occupancy correction is about 5%-7% 
ranging from peripheral to central collisions. The effect of 
occupancy fluctuations within a particular centrality bin 
has also been investigated. The analysis was repeated 
both for a much narrower range in multiplicity in the 
same bin, as well as using the correction weights for 
a higher-multiplicity bin, and the impact on the flnal 
results has been found to be negligible. 



IV. RESULTS 

The flnal 2-D two-particle inclusive correlation func- 
tions for charged particles after all corrections are shown 
in Fig. [5] as a function of A?] and A(/) for the most central 
10% Cu-I-Cu (b) and Au-|-Au (c) coUisions at yfs^ — 
200 GeV. For comparison, the previous measurement in 
p+p collisions (a) at the same energy is also shown. 
In p+p collisions, the complex 2-D correlation structure 
is approximately Gaussian in Ar; and persists over the 
full A(/) range, becoming broader toward larger A(/). This 
feature can be qualitatively described by an independent 
cluster approach, as will be shown later. Compared to 
p+p collisions, heavy ion collisions show not only the 
cluster-like structure, but also a cos(2A(/)) modulation 
due to elliptic flow [l(j|. 

As was done in the analysis of p+p collisions 0, the 
2-D correlation functions in A+A are integrated over 
A0 to give the 1-D correlation functions i?(A7/), shown 
in Fig. [3] This allows a quantitative study of cluster 
properties in pseudorapidity space, with the elliptic flow 
contribution averaging to zero. R(Ar]) is then fltted to 
a functional form derived in Ref. [5| in an independent 
cluster emission model: 



R{At]) = a 



mixed 



(A,/) 



- 1 



(5) 



where p"^^^'"^ (Arj) is the background distribution ob- 
tained by event-mixing and averaged over Ac/). The 
parameter a is equal to Kcs—l, where KcS is the effective 
cluster size, defined as: 



(K) 



'K 



(K) 



(6) 



which depends on the first two moments (mean and 
sigma) of the distribution of cluster size K (number of 
particles decayed from a cluster) over all clusters in all 
events. The function r(Ary) is a normalized Gaussian 
function ^^^ exp(— (A?7)^/ (45^)), where the 6 parameter 
is equal to the width of the Ary distribution of the 
particle pairs from a single cluster. It is connected with 
another variable also characterizing the cluster width, 
CTjj (understood as the width of the distribution of the 

difference T^particic- ■ ^ k,, fi,„ x— IK 



-'7ciustor), by the formula: 5 = \ t^zii 



Gri, for fixed K. The factor 



■j^-j- difference is due to the 

fact that the average of jyparticic in a cluster is constrained 
to be conserved (equal to ?7ciustcr)- Of course, without 
direct knowledge of the distribution of K, the average 




FIG. 3: Two-particle pseudorapidity correlation functions, averaged over the A<j!> range from 0° to 180°, in Cu+Cu (upper row) 
and Au+Au (lower row) collisions for five different centrality classes at y^sj^ = 200 GeV. The solid curves (black) correspond 
to Eq. [5] with the final values of the parameters (see text for discussion). The error bars and bands (grey) correspond to 
point-to-point systematic errors and overall scale errors respectively with 90% C.L. The statistical errors are negligible. 



cluster size (K) and width ct^ cannot be derived based 
on and S. 

Correlation functions for bins in vertex and centrality 
are individually fit using Eq. [5] to extract the effective 
cluster size, Kcs, and the cluster decay width, 5, for 
each bin. These results are then averaged over the vertex 
range to find the final results for each centrality. The 
averaged correlation functions are shown in Fig. [3] along 
with a line showing Eq. [5] with the final averaged values 
of the fit parameters. The three most central points (a 
region of |A?7| < 0.45) in R{Ar]) are excluded from the 
fits mainly due to the large uncertainty stemming from 
residual detector effects. 

Results on effective cluster size (Kcs) and decay width 
(d) as a function of the fractional cross section 1 — cr/ao, 



where ctq is the total A-l-A inelastic cross section, are 
shown in Fig. |4] for Cu-|-Cu and Au-t-Au collisions at 
^s„„ = 200 GeV. The systematic uncertainties are 
estimated using a similar procedure to p-l-p collisions @ 
with an additional contribution from the occupancy 
corrections. The overall scale error, common to both 
Cu-I-Cu and Au-l-Au, is 5% for both K^s and i5. The 
shaded band indicates the value found in y/s = 200 GeV 
p-|-p collisions, which suggests that the cluster properties 
are similar in p-l-p and A-l-A systems. This implies 
that the phenomenological properties of hadronization 
appear to be similar in p-l-p and A-l-A. However, an 
increase of both the effective cluster size and decay 
width is observed going from p-l-p to peripheral A-l-A 
systems. Toward more central collisions, it is also 
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FIG. 4: (a) Kea and (b) 5 as a function of fractional 
cross section {1-a / ao=l for the most central collisions) for 
PHOBOS data (open symbols) and from the AMPT model 
(lines) in Cu+Cu (squares) and Au+Au (circles) collisions for 
|?7| < 3 at ys„„ = 200 GeV. The error bars for data points 
represent systematic errors with 90% C.L. Results from p+p 
collisions at -y/s — 200 GeV 0] are shown by the shaded band. 



observed that the effective cluster size systematically 
decreases with increasing collision centrality in both 
Cu+Cu and Au+Au collisions, whereas the cluster decay 
width is approximately constant over the whole centrality 
range within the systematic uncertainties. Furthermore, 
by comparing the two systems at the same fraction 
of the inelastic cross section (which is related to the 
ratio of impact parameter b to the radius R of the 
nucleus, b/2R), a "geometric scaling" feature is revealed, 
which shows a similar effective cluster size at the same 
collision geometry of the system, i.e. the shape of the 
overlap region. This feature is not obviously expected as 
the cluster parameters are constructed to reflect short- 
range correlations in rapidity and thus are not directly 
connected with the overall geometry of the initial state 
of the collision. Comparison of the data with AMPT [11 1 
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FIG. 5: Near-side (left column) and away-side (right column) 
KcS (upper row) and 5 (lower row) as a function of fractional 
cross section (l-(T/cro=l for the most central collisions) for 
PHOBOS data (open symbols) and from the AMPT model 
(lines) in Cu+Cu (squares) and Au+Au (circles) collisions for 
|?7| < 3 at .ys„„ = 200 GeV. The error bars for data points 
represent systematic errors with 90% C.L. Results from p+p 
collisions at ^/s — 200 GeV @] are shown by the shaded band. 



shows that the model gives the same qualitative trend 
as the data in the same rj acceptance, but with K^s 
values systematically lower by about 0.4. Note that the 
values of K^s and 5 are extracted in a limited acceptance 
of |?7| < 3, and therefore are normally smaller than for 
a full acceptance measurement. The acceptance effect 
will be discussed quantitatively in the next section. In 
AMPT at the generator level, the decrease in effective 
cluster size with increasing event centrality appears to 
be related to the hadronic rescattering stage. Turning 
off hadronic rescattering processes in AMPT leads to a 
larger effective cluster size in both Au+Au and Cu+Cu 
that is approximately invariant for all centralities. 

Further detailed studies on cluster properties have also 
been performed. The entire analysis was repeated for 
pairs in the restricted range. Instead of averaging 
over the whole /S.<j) region, the cluster parameters can 
be extracted for pairs in the near-side and away-side 
ranges (0° < < 90° and 90° < A(/) < 180° 
respectively). Clusters with highp^ generally contribute 
more pairs at the near-side, while away-side pairs mainly 
come from lower px clusters. In this restricted averaging, 
the cos(2A0) elliptic flow component again averages to 
zero. The results are shown in Fig. [S] as a function of 
fractional cross section for Cu+Cu and Au+Au collisions 

y'^NN — 200 GeV. The overall scale error is 3% for 
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the near-side data and 5% for the away-side data, for 
both Kcs and S. Over the studied centrahty range, the 
effective cluster size extracted for pairs from away-side 
decreases by about 30-40% with increasing centrahty, 
whereas the decrease for the near-side is somewhat 
smaUer. Such a behavior could be understood in a 
scenario where the medium is extremely dense at more 
central collisions and only clusters produced close to the 
surface can survive. Then, for low px clusters, it is 
more likely that part of its decay particles travel into 
the medium and get absorbed, resulting in a suppression 
of away-side correlations. As for the observed collision 
geometry scaling of the effective cluster size, it might 
be related to the surface to volume ratio of the system. 
More detailed modeling is still being investigated to 
understand these phenomena. In this case, AMPT shows 
a smaller difference of cluster properties between near 
and away-side than that observed in the data. 



V. CORRECTING FOR LIMITED 
ACCEPTANCE 

As mentioned in the previous section, some particles 
from cluster decay fall outside of the PHOBOS Octagon 
detector acceptance in pseudorapidity {\r]\ < 3). This 
reduces both Kcs and S. In order to quantitatively study 
this effect, a simple independent cluster model (ICM) 
as well as several dynamical models are used. In our 
ICM approach, for each event clusters are generated, 
each with a given mass, transverse momentum, px, and 
ry. The cluster decays isotropically in its rest frame into 
K particles (assumed to be pious) constrained by the 
available phase space. The mass of the clusters is taken 
to be 0.35 X K (GeV/c^). The px and rj of the clusters are 
drawn from distributions that have been tuned such that 
the final inclusive charged particles match the measured 
spectra In addition, the global momentum of 

the clusters is always conserved event- by- event in order 
to preserve the cos(A(/)) component typically seen in the 
p-|-p 2-D correlation function (e.g. in Fig. [5^). 

In Fig. [6^, an example of a 2-D two-particle correlation 
function from ICM is shown with K = 3 measured in an 
acceptance of jryj < 3. It shows a qualitatively similar 
structure to what is observed in p-|-p collisions in Fig. [5^, 
a Gaussian shape short-range correlation along At] which 
gets wider going from near-side to away-side. The narrow 
peak in the near-side part of the correlation function 
is primarily from higher px clusters as well as 3-body 
decays (e.g. in Ref. Q), whereas the broader away-side 
generally arises from the decay of clusters with lower px ■ 
Fig. shows a comparison between 1-D Ary correlation 
functions measured in the limited acceptance of jr;! < 3 



and the full acceptance. From this comparison, one can 
see that the shape is significantly modified by the limited 
acceptance, resulting in reductions in both K^s and S. 
By comparing the cluster parameters extracted for the 
acceptance of \tj\ < 3, if'ff'^^ and jl'^K^, with the full 
acceptance, k]!^^°° and (5l''l<°° an acceptance correction 
can be obtained in the context of the ICM. 

Note that in the ICM, S is not an independent variable, 
but rather depends on the cluster px for given cluster 
mass. The higher the px is, the narrower the cluster 
width will be. In a scenario of isotropic decay of 
the clusters, the maximum possible width turns out to 
correspond to ~ 0.75 for AT = 2 and px of all 

clusters fixed at 0, which is lower than the cluster width 
observed in A-l-A collisions (^ 0.9 in the most peripheral 
collisions). Therefore, to generate a wider range of pos- 
sible 6 parameter values in the ICM, another parameter, 
7, is introduced to artificially modify the width of the 
clusters by rescaling the relative pseudorapidity of each 
decay product: 

= 7 X ('7» = 1,2,3... (7) 

where rjQ represents the pseudorapidity of the original 
clusters, and rji and rj'^ correspond to the pseudorapidity 
of decaying particles from clusters before and after 
modification. In this way, any value of the cluster width 
can be obtained while keeping the original rj of the 
cluster unchanged. For each set of K and 7, the dN/drj 
and dN/dpx distributions of the clusters are tuned to 
match those of the final state inclusive charged particles 
measured in the data. 

The ICM, with a range of values of K and 7, has been 
used to generate a set of two-particle correlation func- 
tions, both for I77I < 3 and the full acceptance. In Fig. [3 
the ratios (A:l'j^l<^-l)/(i^l^l<°°-l) and jl'Jlo/^l'fKoo fj.^^ 
the ICM are shown as a function of 6^"^^^^, as extracted 



(a) Independent Cluster Model (b) 




FIG. 6: (a) 2-D two-particle correlation function in A-q and 
Act) for ICM with A = 3 and 7 = (defined in Eq. g. 
(b) Comparison of 1-D pseudorapidity correlation function 
for ICM with A = 3 and 7 = (defined in Eq. [7} between 
the acceptance of \r]\ < 3 and the full acceptance. 
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FIG. 7: (a) Ratio of K 



\v\<3 
cff 



1 to K™<°° 

cri 



-1 and (b) (jl'^K^' 

to ^I^K"" as a function of i^l^K^ obtained in ICM (solid 
symbols) as well as PYTHIA, HIJING and AMPT models 
(open symbols). The solid line is a smooth function fit to all 
the models. 



directly from fits to the correlation function. It is 
clearly seen that the suppression of these ratios are 
primarily a function of (jl'''^'^ only, and both Kcs and 
S are suppressed more as S^"^^^^ increases. This is 
because correlated particles are more likely to fall outside 
the measured region as the cluster width increases. 
The suppression factors in the dynamical models like 
PYTHIA [15l for p+p, HIJING [3, AMPT for Cu+Cu 
and Au+Au at various centralities are also calculated. 
All models are consistent within about 5%-10%. A 
second order polynomial function is fitted to the values 
of all models in order to generate a smooth correction 

\n\<3 
ff 

and 



function. This function is applied to the measured K^^^ 
and (^'''l^'^ using Eq. [S]in order to estimate K^^^" 
^l»7l<ao £pj. ^j^g experimental data: 



cff 



1) 



MC 



x\v\<oo __ 
data 



r|jj|<oo 

"mc y 

rl')l<3 

"mc 



1) 



MC 



g\v\<3 
data 



data 



(8) 



The scattering of the points around the fitted correc- 
tion curves in Fig.[7]is taken into account as one source of 
the systematic uncertainties on the acceptance correction 
procedure. As a cross check, the cluster parameters in 
the data have also been measured in |?7| < 2. The ratios, 
(if H<2 _ l)/(ifH<3 _ 1) and 5l''l<2/,5l''l<3, are found to 
be consistent with the ICM and dynamical models. The 
residual discrepancies between the results extrapolated 
to full phase space from \r]\ < 2 and \r]\ < 3 are used 
to estimate a separate contribution to the systematic 
uncertainty on the acceptance correction. The total 
uncertainty on the correction is thus found to be 12% 



After applying the acceptance correction, the cluster 
parameters at full phase space, i^),g^°° and (5l''l<°°, in 
p+p, Cu+Cu and Au+Au colhsions at — 200 GeV 

are shown in Fig. [5] for inclusive (left column), near- 
side (middle column) and away-side (right column) along 
with the results from AMPT. The systematic errors 
come not only from the measurement itself but also the 
acceptance correction procedure. The values of KeS and 
6 in p+p collisions extrapolated to full phase space are 
larger than those presented in Ref. @ measured in a 
limited acceptance of jTyj < 3, and better reflect the 
properties of the clusters produced in these reactions. 
Since S measured in both Cu+Cu and Au+Au collisions 
at PHOBOS only weakly depends on centrality, the 
geometric scaling feature of Kcs between the two systems 
still holds after the acceptance correction as shown in 
Fig.m That said, the large values of and S^''^<°° 

clearly pose a challenging question as to the origin of 
such strong correlations with such a long range. In 
calculations from the Therminator model that include 
all known resonances (l7| , Kgg is approximately 2 and 6 
is no larger than 0.75, while in peripheral A+A collisions, 
there appear to be clusters that decay into 5-6 charged 
particles with much larger S. This was not something 
expected from previous data in p+p collisions, although 
there are data on (pr) fluctuations from STAR [ll] and 
PHENIX [T^ that have been interpreted as evidence for 
similarly large clusters in Au+Au collisions 20]. The 
production of jets is a natural mechanism to induce 
clustering phenomena, although one would expect jets to 
lead to a smaller S than isotropic decay. It is also possible 
that additional correlation sources, such as dynamical 
fluctuation of the dN/drj distribution event-by-event, 
may modify the two-particle correlations in a way that 
leads to an increase of the observed effective cluster size. 
The fact that (Sl''l<°° in peripheral A+A collisions is 
larger than p+p, and far exceeds the value expected for 
isotropic decay, also begs the question as to how the 
cluster decays are "elongated" in phase space. Finally, 
it is observed that cluster parameters in central events 
approach values measured in p+p collisions, while those 
in peripheral events are substantially higher — almost a 
factor of two in terms of i^eff — 1 (which is equivalent 
to the so-called "conditional yield" in analyses involving 
high pt triggered hadrons) . Overall, more theoretical in- 
sights are needed to understand these surprising features 
of two-particle correlations in heavy ion collisions. 




FIG. 8: Inclusive (left column), near-side (middle column) and away-side (right column) KcH (upper row) and 5 (lower row) 
as a function of fractional cross section (l-cr/o"o=l for the most central collisions) measured by PHOBOS (open symbols) and 
from the AMPT model (lines) after correcting to the full acceptance in Cu+Cu (squares) and Au+Au (circles) collisions at 
^Sj^j^f = 200 GeV. The error bars for data points represent systematic errors with 90% C.L. Results in p-l-p collisions at 
^/s = 200 GeV after correcting to the full acceptance are shown in the shaded band, which are larger than those in Ref. ^ 
measured in an acceptance of \rj\ < 3. 



VI. SUMMARY AND CONCLUSION 

The two-particle correlation function for inclusive 
charged particles has been extensively studied over a 
broad range in A77 and Acj) in p+p, Cu+Cu and Au-|-Au 
collisions at y^s^ = 200 GeV. In particular, it has 
been shown that the correlation functions in heavy ion 
collisions are not very different from those found in p-l-p, 
allowing a similar interpretation in terms of clusters. In 
this approach, multiple particles are understood to be 
emitted close together in phase space, with a typical 
effective cluster size of 2.5-3.5 charged particles in p-fp 
collisions. The correlation functions in A+A show a non- 
trivial decrease in effective cluster size with increasing 
centrality, and a surprising geometric scaling between 
Cu-|-Cu and Au-f-Au collisions. Analysis of near- and 
away-side clusters provides additional information on 



the details of the cluster properties. Extrapolating the 
measured cluster parameters to the full phase space using 
an independent cluster model as well as other dynamical 
models such as PYTHIA, HIJING and AMPT, the 
effective cluster size and width increase in magnitude 
to a level which seems to challenge most conventional 
scenarios of the hadronization process. Clearly, more 
experimental and theoretical work will be needed to 
understand these novel aspects of heavy ion collisions. 

This work was partially supported by U.S. DOE 
grants DE-AC02-98CH10886, DE-FG02-93ER40802, 
DE-FG02-94ER40818, DE-FG02-94ER40865, DE-FG02- 
99ER41099, and DE-AC02-06CH11357, by U.S. NSF 
grants 9603486, 0072204, and 0245011, by Polish 
MNiSW grant N N202 282234 (2008-2010), by NSC 
of Taiwan Contract NSC 89-2112-M-008-024, and by 
Hungarian OTKA grant (F 049823). 



[1] R. E. Ansorge et ai, UA5 Collaboration, Z. Phys. 
Particle and Fields C 37, 191 (1988). 



[2] F. Henyey, Phys. Lett. B 45, 469 (1973), 
E. L. Berger, Nucl. Phys. B 85, 61 (1975), 



10 



J. L. Meunier, G. Plaut, Nucl. Phys. B 87, 74 (1975), 

C. Michael, Nucl. Phys. B 103, 296 (1976). 
[3] K. Eggert et al, Nucl. Phys. B 86, 201 (1975). 
[4] D. Drijard et al, Nucl. Phys. B 155, 269 (1979). 
[5] A. Morel and G. Plaut, Nucl. Phys. B 78, 541 (1974). 
[6] B. Alver et al, Phys. Rev. C 75, 054913 (2007). 
[7] L. J. Shi and S. Jeon, Phys. Rev. C 72, 034904 (2003). 
[8] B. B. Back et al, Nucl. Inst. Meth. A 499, 603 (2003). 
[9] B. B. Back et al, Phys. Rev. Lett. 87, 102303 (2001), 

Pradeep Sarin, PhD thesis, Massachusetts Institute of 

Technology (2003). 

Wei Li, PhD thesis, Massachusetts Institute of Technol- 
ogy (2009). 

[10] N. Borghini, P. M. Dinh, and J.Y. OUitrauh, Phys. Rev. 

C 63, 054906 (2001). 
[11] Z. W. Lin, G. M. Ko, B. A. Li, B. Zhang, and S. Pal, 

Phys. Rev. C 72, 064901 (2005). 
[12] B. B. Back et al, Phys. Rev. Lett. 91, 052303 (2003). 



[13] B. Alver et aZ.. larXiv :0709.4008yl (2007). 

[14] S. Adler et al, Phys. Rev. C 69, 034909 (2004). 

[15] T. Sjostrand, S. Mrenna, P. Skands, JHEP. 0605, 026 

(2006). Version 6.325, single diffractive process excluded. 
[16] M. Gyulassy and X. N. Wang, Gomput. Phys. Gommun. 

B 83, 307 (1994). Version 1.383, single diffractive process 

excluded. 

[17] A. Kisiel, T. Taluc, W. Broniowski and W. 

Florkowski, Gomput. Phys. Gommun. 174, 669 (2006) 
arXivmucl-th/0504047 . Calculation by A. Kisiel, 

private communication. 
[18] J. Adams et al, Phys. Rev. C 71, 064906 (2005), 

J. Adams et al, Phys. Rev. C 72, 044902 (2005). 
[19] K. Adcox et al, Phys. Rev. C 66, 024901 (2002), 

S. S. Adler et al, Phys. Rev. Lett. 93 092301 (2004). 
[20] W. Broniowski, P. Bozek, W. Florkowski and B. Hiller, 

PoS CFRNC2006, 020 (2006) arXivmucl-th/061 1069| . 



